###########################################################################################
# 
# Prepare List Experiments and Test for No Design Effects
#
########################################################################################

### Democratic Republic of Congo

table(congo$d_listCRSV)
table(congo$list1_CRSV)

congo$treatment <- congo$d_listCRSV
congo$Y <- congo$list1_CRSV

ind.congo <-lm(Y ~ treatment, data=congo) # List Experiment
c.ind.congo <- coef(ind.congo)[2]
se.ind.congo <- se.coef(ind.congo)[2]

design.congo <- ict.test(congo$Y, congo$treatment, J=3) # Test for No Design Effect

########################################################################################

### Liberia

table(liberia$cw_random)
table(liberia$cw_list_count_control) 
table(liberia$cw_list_count_treat)

liberia$treatment <- ifelse(liberia$cw_random>.5, 1, 0)
liberia$Y <- ifelse(liberia$treatment==1, liberia$cw_list_count_treat, liberia$cw_list_count_control)

ind.liberia <-lm(Y ~ treatment, data=liberia) # List Experiment
c.ind.liberia <- coef(ind.liberia)[2]
se.ind.liberia <-se.coef(ind.liberia)[2]

design.liberia <- ict.test(liberia$Y[is.na(liberia$Y)==F], liberia$treatment[is.na(liberia$Y)==F], J=3) # Test for No Design Effect

########################################################################################

### Sri Lanka

sri$treatment <- ifelse(sri$Questionnaire==2, 1, 0)

sri$Y <- ifelse(sri$Questionnaire==2, sri$J7, NA) 
sri$Y <- ifelse(sri$Questionnaire==1, sri$I7, sri$J7)

ind.sri <-lm(Y ~ treatment, data=sri) # List Experiment
c.ind.sri <- coef(ind.sri)[2]
se.ind.sri <- se.coef(ind.sri)[2]

design.sri <- ict.test(sri$Y, sri$treatment, J=3) # Test for No Design Effect

########################################################################################
